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1. Introduction 



, In the recent years a considerable progress has been made in investigation of intcgrable models of two- 

' dimensional quantum field theory. S matrices^ and form factors^ have been found in the framework of 

, the bootstrap approach. A free field representation'^ for the form factors of the sine-Gordon and Thirring 

models was proposed. But these achievements are based on some guesses and we need some additional 
physical grounds to understand better the nature of these results. The main problem can be formulated in 
pi • the following way. It is possible to construct sets of form factors, that describe some local operators on the 

D , basis of the axiomatic approach. Moreover, it is possible construct the whole space of local operators. But 

i-^ ' it turns out very difficult to understand which operator corresponds to a given set of the form factors, and 

vice versa: if we consider a given operator, then which set of axiomatically allowed form factors corresponds 
KJi ; to it? 

1*^ , A approach to lattice models of statistical mechanics based on the free field representation of correlation 

■ functions and form factors has been developed by the Kyoto group. '''^'^ In contrast to the situation in the 

quantum field theory this approach has a firm physical basis. There is a simple rule that allows one to 
construct a bosonic representative to any local lattice object. So it one can easily construct a set of form 
factors corresponding to each local lattice insertion. On the other hand, the lattice models are known to be 
described by models of quantum field theory at large scale in the vicinity of the critical point. So, it seems 
promising to investigate these models in the vicinity of the critical point to extract some information on the 
respective quantum field models. Our main strategy is the following. We consider a lattice theory and find 
form factors of some local insertions in the scaling limit. Then from some additional information, such as 
scaling dimensions, conservation laws, symmetries, mutual locality etc. we associate a local operator in the 
field theory to the insertion. This allows us to find to relate form factors to particular operators. Here we 
do not fulfil this program sequentially, but make some steps in this direction. 

In this letter we consider the six vertex model'' in the antiferroelectric region in the vicinity of the critical 
point. In the critical region this model is described by the SU(2) invariant Thirring model or, equivalently, 
by the sine-Gordon model at the critical coupling.^ We analyze the spectrum, the S matrix, and the form 
factors of the polarization operator . Our result does not match that obtained by Jimbo and Miwa, ^ but 
we shall see that it is consistent with the formula for form factors in the sine- Gordon model, guessed recently 
by Lukyanov. ® We show the crucial role of the zero mode operators to the behavior of the form factors in 
the scaling limit, and show how to pass to a 'conjugate' pair of the zero mode operators simplifies the limit. 
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2. Spectrum and S Matrix 

The six vertex model in the antiferroelectric region is described by two parameters x and ( , 

0<x<l, l<C<x~'^, 
so that the weight matrix of the model is given by the trigonometric R matrix 

R{Q+.-R{Q.+ -^ (C) , 4C)-C (^4^_i.^4)^(^2<..^4)^ 



(2.1) 



where {z;p)oo = n^o(l ~ ^p"). The i? matrix satisfies the Yang-Baxter equation, and ( serves as a 
multiplicative spectral parameter. In particular, the row-to-row transfer matrices T{Q = T{x; Q commute 
for different values of C but fixed x. 

The eigenvectors of the transfer matrix at the infinite lattice are 'A'' -particle' states |£i^i, . . . , snZn) , 

N 

T{C)\£izi, . . .,£nzn) = W T{zj/C,)\eizi, . . .,snzn) (2.2) 

i=i 

with £j = ± = ±1 being the 'charge' of 'jth particle', and Zj is the 'rapidity' of 'j th particle', 

\zj\ = 1. (2.3) 
The 'one-particle' eigenvalue function t{z) is given by 



(2.4) 



where Qp{z) — {z]p)oo{p/ z]p)oo{p\p)oo is a multiplicative form of the 9i function with the Jacobi parameter 

Qp{z) = -2ip-^'^z^'^ei{v;T), z = B^^'\ e'"^=pi/2. 



Note that this form of the eigenvalue function is universal for all lattice models based on the A\ type R 
matrices. 

We know that in the scaling limit, when we are bringing the system closer to the critical point, simultane- 
ously enlarging the scale, we have to obtain a system of relativistic particles with the states |£i^i, . . . , £nSn) , 
the Hamiltonian H , and momentum operator P , such that 

AT 



P \siOi, .. ., £nOn) = M^'^shOj \eiOi, .. ., snOn), 



'7 (2-5) 
H\ei6i,...,eN0N) = M^chOj \ei6i, . . . ,£nOn), 

with M being the mass of the particles, and 9j being rapidities. Here the parameter is connected 
somehow with z, and the functions M ch.9 and Msh0 must be related somehow with the function t{z). 
The scattering matrices of the lattice models must as well become the relativistic scattering matrices in the 
scaling limit. 

The scaling limit corresponds to the limit a; ^ 1 , so we set 

x = e~^, e — > 0. 
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Let us introduce also additive spectral parameters 7 and 1? that survive in the scaling limit: 



7,i?eR, 0<7<7r. 



(2.6) 



It is convenient to pass to the conjugate module in the theta functions. Namely, for the function T{'d) 
we have 

The Jacobi parameter, corresponding to the theta functions is 

q = e-"'/2^ . 

We also use the usual designations: k and k' are module and conjugate module, K and K' are half-periods. 

The Jacobi parameter q tends to zero in the scaling limit and wc can use the standard expansions in q 
for theta functions. But before doing it let us understand, which values of the spectral parameter correspond 
to the low-lying excitations which survive at a large scale. The low-lying excitation correspond to the largest 
in absolute value transfer matrix eigenvalues. The physical values of the spectral parameters C and z satisfy 
the constraints I < C < and jz] = 1. So we ought to find the maximum of \t{z/Q\'^ in z on the unit 
circle for fixed C in the interval 1 < C , or, in other words, the max;imum of |t('!? -|- i7)|^ in real for 

fixed 7, < 7 < TT. To do it, note that the function 

\l \" -r — , ^ 1. i7r\ /I / ii? , 7 I 1. i7r\ ^ ' > 

^4 l2¥ 2^ 4' 2ij ^4 (^2^ 2^ + 4' 2ij 

of 1? can be analytically continued to a doubly periodic function on the complex plane. Then its derivative in 
'd has two double poles in the periodicity region, and, therefore, has exactly 4 zeros. Using the symmetries 
of this function it is easy to check that all 4 extremal values are situated at the 'symmetric points' and two 
of them on the real axis. The maximum is achieved at the point i? = and the minimum at i? = 7r^/2e. 

We see that we ought to consider the limit of the elliptic functions q ^ Q (K— !-7r/2, K'— >oo) for the 
values of i? in any finite interval around 0. Applying to (2.7) the standard limiting formulas 

dn = [i+4gcos2a + 0((z2)], 

TT Z±V 

fc' = l_8g + 0(g2), K = |[l + 4g + 0(g2)], 

we obtain 

r(i9) = l+4igshi? + 0(g2). 
The system can be considered continuous if 

Let us think that our system lives on the lattice {am,bn), m,n G Z, a and b being lattice parameters. 
We will think of am as of the imaginary time it' and of bn as of the space coordinate x' . Recall that the 
operator T(l) is the shift operator by 6.^ So the momentum operator p' is given by 

eip'Wb ~ 1-H4igshi? 

or 

p'(^) = ^shi?. 
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The Hamiltonian is given by e = T{z) , and for the one-particle energy we have 

^-E'{^)a ^ + i^) ^ 1 + 4igsh('!? + i7) 

or 

4ig Aq 
= sh(t? + i7) = — (sin7ch'!? — icos7shi?). 

To restore the relativistic spectrum we recall that the initial lattice system is essentially anisotropic. Hence, 

we need to change coordinates (t', x') (t, x) , so that the spectrum would take the standard form E{6) = 
M ch.9 , p{9) = MshO . It is achieved by the coordinate transformation 



t = t'^sin7 =— i6sin7-m, 

X = it'^ cos 7 + x' = +b cos j ■ m + bn, 



' b „„„ ^, , ^1 , (2-8) 



and by setting 

e = 'd, ^ = y- (2-9) 

Consider now the S matrix. For the six vertex model it is given by^ 

S{z) = -R{z-^). 

In the scaling limit we obtain 

s{e)XX = sie)zz = -Soie), 

s{9)Xz = s{e)zl = -s,{e)-. ^ 



ITT 



r(i-f)r(i + f) 

This S matrix coincides with the critical sine-Gordon S matrix." 
The sine-Gordon model is the model with the action 



/ 



The dimensionality of the constant A is M^~^^^ . For < /3 < 1 the term cos P^p is a relevant perturbation, 
and the theory is superrenormalizable. For /3 > 1 the perturbation is irrelevant, and the properties of the 
theory are completely unknown. The critical sine-Gordon theory corresponds to the marginal value /3 = 1 . 
In fact, dimensional transmutation takes place at this point, and the theory is an integrable theory of massive 
interactive charged solitons. S{9) is the pair soliton-soliton S matrix. 



3. Free Field Representation £ind Conjugate Zero Mode 

Now let us turn to the calculation of form factors. We start from the free field representation of Ref. 6. 
Consider the Heisenberg algebra 

ro^i ■ r 1 . [2m]^[m]. 



m 

The spaces Hi, i = 0,1 are generated by the operators a_„, n > 0, and exp(±i\/2(5) from the vacuum 
vectors |0)(*^ , defined as 

a„|0)W =0, n > 0, 
^|0)« = ^|0)«. 



" In the case of the six vertex model we do not need the duality transformation here. 
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We also shall use the notation \p) for the vector 

|p)=e^f«|0)W, P\p)=p\p). 
The vertex operators (C) > Vs^~^'^\z): Hi — » are defined as^ ^ 

F|'-''')(C) = Gi/V/4:e'^(«:, 

Jc, 27ri 



y(l-M)(^)^C-l/2^3/4^1/2.gi0(.).^ 



(3.1) 



where 



and 



.a;-|"l/2_2-"_ 



(3.3) 



The normal ordering operation :. . .: places P to the right of Q, and a„ to the right of a_„, n > 0. The 
contours Ci and C2 encircle the point so that x( is inside and x~^( is outside Ci, and is inside 

and xz is outside C2 ■ 

The form factors of local operators O are given by® 

(0|O|£iCi..-^ivCiv)!^ = -1-Tr„. (x4^'^'6(^-^0t;(^a-O(^^)...V;(i-M)(^^)) , (3.4) 

Ail'*' J ^ ' 

Here 

n 



is the shift operator 



and Xi(i) is the character of Hi 



X.(i) = TV„,i^'"=(iV2.t)-i. 



The operator O^''*'^ is determined by the operator O and consists of the operators V^^ *''HC)- We shall 
co^ider here an important example of the polarisation operator = ( J _« 

(hi) 

vertical links of the lattice. The respective operator (cr^) is given by 



b In terms of Ref. 6 yii"'-*) (C) ~ (C^/^) , V'i'"*'*'^ ~ 



Consider an arbitrary trace 



where 



= e-'Oe'^'exp^^; yl„«_„^ '^p(Y. ■*-"«'.)■ 



It is easy to check that the trace can be factorized, 



^(0 _ rrii) 



-njY -Tjv, 



into the contribution of the zero mode {Q, P) , H^^ , and the contribution of the oscillators a„ , Fn ■ The 
oscillators contribution is given by 



Fn = 



logCj 



3 3<3' 
1 



logi^j 



1 



(3.6) 



n' 



Tr* (a; 



4D. 



Here asterisks mean that the contribution of oscillators is only taken into account; (^+ and (j)- are contri- 
butions into (j) from a_„ and a„ respectively (n > 0). 

The zero mode contribution is the most interesting part. Let us write it down explicitly. It is only 
nonzero for J2j '^j = • Let us introduce the notations 



N 



N 



(3.7) 



Then 



3<3' 



neZ+i/2 



n 



V2 



= e---^ E exp(-4e(n^-|)+enV2$:/3,)=e-'^^^3-.(-ie^:^;if) 



i6Z+i/2 

ijfe-'^^exp 



4e 



1 B 



4 4^/2 



73+1 



1 B .71 



4 4V2' 4ey ' 



(3.8) 



where the subscript means the contribution of the zero mode. We passed to a conjugate module in the 
last equality. The most interesting feature of this expression is that we can rewrite the last theta function 
as a trace of an operator. Expanding it in the series, we can easily check that 



H}^Xo{x ) = - W-exp 



2 V e 

xE 

nez 



4 4V2J 



(3.9) 



Here P' and Q' is a new pair of canonically conjugate operators and \p)' designates the eigenvector of P' : 



[P',Q'] 



P'\P)'=P\P)', b)' = e'f« |0)'. 



(3.10) 



Naively, the tcaiii with n = only survives in the limit e ^ . But as we shall see below in the calculation 
of the form factors of this term vanishes, so that the terms with n = ±1 become leading. 

Note also that we can substitute P — > iP'/e and Q — > —ieQ' directly in the vertex operators. This 
substitution does not affect the commutation relations of the vertex operators, but the summation rule over 
the momentum values changes drastically and the operator D^*' must be substituted by 



2^2 



OO 



(3.11) 



The explicit form of the nonzero 'two-point' traces is 
1 



-logFii(7i -72) 

-l0gF22(^l -^2) 

-logFi2(7i - 1612) 



^Tr, {x^'''MCi)MC2)) = J2 ^^^ch2en ( 1 



n=l 



n sh 2en 



^ e"' shen , ^ , .(7i 
> = ch 2en 1 — 1- 



^Tr.(.--^.(.,)^*(.2)) 
^TV,(.-^.(C0^.(.2))=E^"'^"" 



7i - 72 



^ n sh-'2en 

n— 1 



n 



sh^ 2en 



ch2en 1 



71 - 1^2 



Similarly, the traces of the form x* ^{x"^) Tr*(a;'*^*0+0_) are 



(3.12) 



-log[x-'/^G-y^C+] = J2 



■log[ar-=^/4G-V2Ffi(t;,C)C_(7,M)] 



n=l 

00 



2nsh 2en 



-Sen 



sh 



en 



^ 2nsh^2 

n— 1 



-log[x-3/^GV2C+] = E 

00 



en 

e"*^" sh en 
2n sh^ 2en ' 



-4 (^ch 2en^ - l) ch en + (2 ch en - 1^ 



(3.13) 



'sh( 



n=l 



2nsh 2en 
V = e-' ^' '"''^ , w -- 
It is convenient also to introduce the functions 

Ai(7) = Fr.Hl - 7r/2)Ff/(7 + 7r/2), 
^2(7) = f'r2'(7 - 7r/2)Ff/(7 + 7r/2), 
F22(^) = F^-^\e + m/2)F^^\e - m/2), 



-4 ^ch2en 1 ) chen + (2chen - 1)^ 



^11(7) = A"i'(7 - ^2)^1' (7 + 7r/2), 
Fi2(7) = A"2'(7 - 7^/2)^1-21(7 + 7r/2), 
#22 (^) = F2-2i(^ + i7r/2)F2-2i(^ - i7r/2). 



(3.14) 



All traces of the vertex operators can be written in terms of the functions (3.12-3.14). To do it we need 

1 /2 ~ ~ 

to write a function , for each V^iQ) , Ve{zj), a function Fii(7j — 7j') for each pair V+{Q) and 

V+iCf), 3 > /, a function i^ii(7j - 7j')^ii(7j - 7/ - Mj') for each pair ^+(0) and F_(Cj')> J > ete., 
then to write the overall coefficient from the zero mode, and finally to integrate over all variables /Xj and Vj . 



4. Scaling Limit of Form Factors 

Consider the scaling limit of the functions (3.12). Applying the formula 



n=l 



^ \ en^ 



a I e 



+ ef{en)j = —02 - Oi log e + J dij /(r/) + 0{e) 



for /(jj) being finite as — > 0, we easily obtain 



Fii(7) = e-'/V/''(/ii(7)+0(e)), 
F22W=eVV/^'(/22W+0(e)), 



-^12(7) = 9 



— 



exp ■ 



7 
12 



+ 0(e°°) 



(4.1a) 
(4.1b) 

(4.1c) 



with 



■log/ii(7) = 

-log /22(e) = / 
Jo 



drj 

V 
drj 
rj 



sh^ 2r] 

'^^cli2r 
sh^ 27] 



-) 


1 






4:1] 






1 















(4.2) 



The expression for F12 demands a more detailed consideration. It is easy to check, that the formula is 

correct up to all orders in e , because in the decomposition of the summand in e the odd powers of n vanish. 
And it is easy to check that X]^i(— 1)""^?^^^ = 0, p > 1. However, we need a more strong (and more 
rigorous) estimation than 0{e°°). We can give it for the ratio 



^12(7 + ^2) Qx<-xy) 



A2(7-V2) Qx^{-x^y) 



y^'^[l + 0{q)], y = x 



2(1-7A) 



(4.3) 



The last equality follows from the earlier consideration of the function t{z) . We shall see that this estimation 
is strong enough to find the form factors of . 
For C± , C± we find 



^;C_(7,m) ^ e-VV/'^c_^/i(/z)e^T/-, 
C+ 



(4.4) 



with 



logc+ 
log c_ 
logc+ 
log c_ 
/i(m) 

/2M 



drj 
V 

drj 
1] 

drj 
V 

dr] 



shry 
sh^ 2r] 



1 1 

4:1] 4 



e 3^ jJZL(2chr/- 



2 ch 7] 
1 1 



sh^ 27] ' 
sht] 
sh^ 2rj 



1)- 77 + 76-" + 



47] 



2chr? 



4?7 4 



2 Jo V 



e-vJ^(^2chr]-lf 



2chr? 
1 1 



sh^ 2r/ ' 



41] 



exp 



(/z-7r/2)(/x + 7r/2) 

TT 

~(z^ + i7r/2)(i^-i7r/2) 



/ — -tt;- (ch2ry- - 1) 
Jo ?? sh2r? V TT / 



exp / — 



sh2r? 

dr] 



rj sh 2rj 



w 

ch2?7 — 

TT 



2ch77 



1 



(4.5) 



All equalities in (4.4) are up to 0(e) , but the expression for C_ includes the precise dependence on 7, which 
will be important later. 

Wc need to make one more remark. The exact formulas for form factors on the lattice contain integrations 
over closed paths encircling zero. It means that we must integrate over ji and v along to infinite contours. 
But at \v\ ~ 1/e the formulas (4.1-4.5) is no more valid. We should analyze its contribution. We can 
argue in the following way. It is easy to check that 



/ii(ie) 
/22(e) 



|-V4g|e|/4^ 
|V4g|e|/4 



HO) 



\e\- 

|^|-3/2 



1/2 e-l^l/2. 



-\e\/2 
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for large real 9. Therefore for large imaginary fi or large real z/ a factor ^ ^ e^'^' or 

~ f22 {^) ~ |z^|~^e~l'^l appears. We shall see that the zero mode gives the factors ~ e'^' , e''^! in the 
case of (j^ . So the contribution of i/j 1/e is negligible, but the contribution of 1/e must be taken into 
account. Nevertheless, as we shall discuss later, the correct answer can be extracted from the approximation 
for small /x. 

At last, we define the states 

|£l6'i, . . . ,£jv6'jv)|n = ( — ) kl^l, • • ■,SNZN)ln 



so that it would have the natural for two-dimensional field theory normalization: 

TVT / (27r)^ ini^l^l' ■ ■ ■ ,£N0'N\ei9i, . . . ,£jvgjv)|n = 1- 

Now we are ready to calculate of the form factors of in the scaling limit. The form factors of 
read precisely 

{Q\a^\e,e„...,e2ne2n)[:^ = "X' ^X' (^F ^ "^-L(7,/^;£,^,^)^^-aa(7,/^;£,^,^^), (4.6) 

where H and F are contributions of the zero mode and the oscillators respectively. Let us begin with the 
oscillators. 

The contour C[ goes from — ioo to +ioo to the right of —^tt and to the left of ^tt, and the contour 
C2 goes from —00 to +00 above ^tt and below —^w (with an inflection). 
The oscillators contribution factorizes as follows: 

2e 

F-oca{l,tJ.;e,0,v) = —F_aa{l^)F{£,6,v)S-aa{l,lJ-;£-,0-,v). 
TT 

Here F_aa (m) is the F function for the vacuum expectation of : 



(4.7) 



J^-+(m) = -C_(7 + 7r,M)C+Fn(7r)i?n(/i + 7r), 

TT 

F+_(m) = -C+C_(7,M)^^ii('r)Ai(7r-M). 

TT 

^From the evident property ^11(7 + 2'7r) = -Fii(— 7) we obtain 

^^_+(M)-C_(7 + 7r) • ^^-'^ 

Note that this result is precise. Up to 0{e) the quantities _F_+(^) and F+-{n) coincide and 

F_+ip) ~ c+c_/i(Ai)/n(^)/n(/i + n). (4.9) 

The function F{£, 6, v) is the factor that only depends on the asymptotic states rather than on the local 
operator: 

2n 



j6JV_ j<j' 

n F22{9r-9,+i^,.) n F22{9r - 9, - n F22{9r-9, + i^,.-i^,) 

-j n /2(^.>^.)n/22(^/-^i) 



X n f22{0j'-9j + iyj,) n /22(%-^i-'-i) 

X n /22(^,'-^,+I^/-I^,) (4.10) 



with N± = {j\sj = ±} . From the requirement aj = we have Sj = or #-A^+ = #N- . 
The third factor S'_Q,a(7, /x; e, 9, v) is given by 

S-+ (7, m; £, ^, = n (7 - i^j + T^)Fi2 (7 - i^) A2 (7 - i^j + M + tt) 

i 

X -^12 (7 - i^'j - ii'j + 7i')Fi2(7 - i6'j - i^'j)-F'i2(7 - i^'j + M - '^^0 + ti"), 

5+_ (7, m; £, ^, i^) = n ^12 (7 - i^i + 7r)Fi2 (7 - i^) A2 (7 - i^i + m) 

i 

X W ^12(7 - i^j - ii'j + 7i")A2(7 - i^j - i^'j)-F'i2(7 - ^^j + M - i^^j), 
and it is equal to 1 up to 0(e) . But the ratio can be estimated much more precisely by use of Eq. (4.3): 



(4.11) 



n 



n 



Fi2(7 - \9j + M - i^'j] 



^ F12 (7 - 16 J + /X + tt) F12 (7 - i% + /X - iz^j- + tt) 



ei^®A + 0(g), 



(4.12) 



Note that 8 = 'mB' /\/2 with B' being the contribution of into B from Eq. (3.7). 

We conclude that the factor F-aail, l^',^,^,^) is independent of a up to 0(e), but we are able to 
calculate the ratio F+-/F-+ with much better precision, namely up to 0{q). This precision allows us to 
treat the terms with n = ±1 in Eq. (3.9), which are of the order g^/^ 

Now consider the zero mode contribution H-aa ■ From Eq. (3.8) we obtain 



^fiL(7,M;^,^,-)=e--'-(^-)e-^^ 



exp 



4e 



1 - a 



1-a B' 
—— + 



4\/2 27r'4e 



TT"! :r ) / ( 7; :r ) ' 



B' n 



1 Itt 



4' 4e 



(4.13) 



where A' and B' are contributions of 14 into the constants A and B from Eq. (3.7). Collecting equations 
(4.8), (4.12), and (4.13) we obtain for the functions in the integrand of Eq. (4.6) 



H+-F+. 



f B' ij, m 
\AV2 ~ 2^' ii 



/ B' IJ, in 
\AV2 ~ 2^' ii 



+ 0{q) 



(4.14) 



[recall that 94{u) = 93{±l/2 + u)]. It means that the difference F_+ — f+_ is proportional to the difference 
of theta functions in the scaling limit. 



In our case 



6'3+i(u;i7r/4e) - e4_«(w;i7r/4e) = (-)*g^/^ • 4cos27rM + 0(9). 
B' 



(4.15) 



27r aV2 



1 
'2^ 



We see that the contributions of the order g'^ [terms with n = in Eq. (3.9)] cancel each other in cr^ , and all 
form factors are comparable and proportional to q^^^ ~ M^/^ _ Once we extracted the factor of this order, 
we may only take into account the terms of the order e*^ in all other multipliers. We obtain the expression 



(O|a^(O)|ei0i,...,£2„M 



)cos( /X- -6 ), 



(4.16) 
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where F !-(//,) and F{e, 0, v) arc given by approximate equations (4.9) and (4.10). Wc sec that the integrand 

nearly splits into two factors, one of which only depends on fx and the other one on uj . The only factor 
depending on both and vj is the 'zero mode' factor cos(/i— ^O). Note that the function F_+{/j,) is even. 
Therefore, wc may substitute cos(/U — ^0) in the integrand by cos/i cos^©, and the integral factorizes 
completely. Consider the function F__|_(/i) cos/i. Calculating it numerically by use of the Mathematica 
package I found it to be a constant: 

F_+(/i) cos/i = —1. 

Hence the integral in /x is linearly divergent. This is the problem mentioned above. Therefore we should 
consider large e~^. Note, that Eqs. (4.1c), (4.3), (4.8), (4.12), and therefore (4.14) are valid in the 

case of large /U. The equation (4.15) is valid in the whole region — 7r^/2e < Im/x < 7r^/2e except small 
vicinities of the ends. To find approximations of Fn , C_ in the region ;» 1 one can expand summands 
in en everywhere except the places where en enters with a large coefficient n, then one can sum the 
series in the regions of convergence and continue analytically the answers. Such consideration shows that 
the integrands in all form factors of remain factorable and constant in fi in nearly the whole region 
— 7r^/2e < Im/x < 7r^/2e. On the other hand the integration region in n for finite e must be precisely 
[-i7rV2e,i7rV2e]. 

So we finally obtain 

(0|a^(0)|Mi,...,£2n^2„)[^^ 

= (_)^2 / 4 \ r jp^p^ ^ ^^^^ 

^ ^ TT^ ' ^ \bM J Jc'^ (27r)» ^ ' ' ^ 2 

- <-'■ iT {^^ k (S- <^ - ^ I,}'' ' "'^') 

X n /2(l^) n /22(^,' - ^i) n /22(^,' - + Vr) n /22(^,' - - V^) 

jeN- j<j' j<j',j'eN- j<j'jeN- 

X n /22(^/-^,+^/-^,). (4.17) 

For n = it matches the well known Baxter-Kelland formula 

Let us interpret the result from the point of view of quantum field theory. First of all, note that every 

form factor is proportional to M^^"^ . It means that the scaling dimension of the field a^{x) is 1/2. There 
are two linearly independent local fields of this dimension in the critical sine-Gordon theory, e"''/^ and 
e-'v/2 ■\Yg expect that the field cr^{x) is a linear combination of these fields. Comparison with the result 
of Ref. 9 supports this conclusion and makes it possible to calculate coefficients of the combination. The 
decomposition of corresponds to the decomposition of chi0 into two exponentials in Eq. (4.17). To 
make it clear, let us pass from the bosonization of the lattice model to the bosonization of the field theory. 

The 'conjugate zero mode' operators P' and Q' survive in the scaling limit, but the coefficient at P'^ 
in tends to infinity. So we may only take into account P' . The only trace of Q' lies in the requirement 
J2j 0!j = . Considering form factors of lattice objects like we must accurately extract the leading 
terms in (3.9), as we have done above. But for the field theory operators we can consider any expectation 
(P'l . . . . For e ^ it takes the form 

if«(P') = (6M/4)4^"/-e'(-i)'^^' e'/^"^' . . . e"^^^' . (4.18) 

It means that the zero mode contribution to the scaling dimension of the fields is equal to 4P'^/7r^. The 
contribution (4.18) just provides the factors before exponent field in Lukyanov's formula for vertex operators 
corresponding to the field e""^ for P' = 7ra/V2. Note that in this case AP''^/-k'^ = 2a? gives the correct 
conformal dimension of the field e""'' . 
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Consider the oscillators a„ . It is tempting to think that the set of a„ tends in the scaling limit to a 
continuous set of operators a{r]). But it is not quite correct. To see it, let us consider the field ^12(7). The 
expression for this correlation function contains the factor (—1)" [see Eq. (3.12)], which is rapidly oscillating 
in the scaling limit for finite 7 . It reflects the fact that the oscillators forming V± differ from those forming 
V± by the factor (—1)" [see Eqs. (3.1), (3.3)]. It means that the set of a„ splits in two continuous sets of 
operators a{ri) and b{i]), 

lain), a{v')] = Mv), b{v')] = 5{v + rj'fJ^^^, [a(r?), %')] = 0. (4.19) 

Now let us define the operators 4>{j) and (t>(9) as follows'^ 

(j)M = y^P'j - f dr] — ^ e-l"!/^ e 

TT J_ 



V2„,., r .^ /^fa) ,-1.1/2 

ish 27) 

(4.20) 



Finally the definition of the vertex operators is the following 

y+(7)=.9i/2:e''^(-'):, 



1/2 / ^ ^ . i0(7)-i./'(7+M-7r/2)-i,ji(7+M+7r/2). 



V (9) - -,/l„-i/2 f ^ ^ i^(g)-ij(e+^+i^/2)-i0(e+i.-i^/2). 



with 



log 5 



= 

Jo V 



dr] e ^ {chij — 1) 



2chr? 



if we assume the regularization rule 



TT 

const. = — const. 
2e 



for f{r]) finite at 77 — > 0. 

The vertex operators V±{9) are Lukyanov's vertex operators from Ref. 9. The factors before the 
exponential operators are hidden in the zero mode P' . The operators V±{'-f) can be considered as the 
primed operators from Ref. 3 at the infinitely distant point. They only interact with the operators V±{6) 
through the zero mode. Their physical meaning in the framework of the quantum field theory is vague. As 
we can see, the initial idea by Lukyanov that the primed vertex operators describe local fields was in the 
right direction, but it is necessary to move them to the point 7r'^/2e which is infinitely far for e — > 0. 

Comparison with Ref. 9 gives 

cr^(.T) - cos i</7(.x). (4.22) 

Unfortunately we are unable to find the proportionality coefficient from the vacuum expectations, found 
by Lukyanov in Ref. 11. In the limit /3 — > 1 these expectation values become infinite. It means that the 
critical sine-Cordon is not described by a free boson in the ultraviolet region, and the normalization of 
the operators of Ref. 11, based on the conformal field theory, is not applicable in this case. In the lattice 
theory it is reflected by the additional factor ~ log(4/6M) in the form factors. This factor means that 



The terms ~ logo;'' in (3.3) do not contribute traces in this hmit. In the exact formulas they can be also omitted together 
with the factors x^/^ in (3.1). 
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the normahzation changes with the scale. Maybe it witnesses about logarithmic factors in the correlation 
functions at short distances. 

Our last remark concerns the fields^'^'^ cr* = ^{'^^ ^ iu^). The leading zero mode contribution to 
these fields comes from n = in (3.9), but the scaling dimension 1/2 comes from the oscillators, because 
the requirement J^j '^j — gives here — = ±2, and the factors depending on q in (4.1) and 

(4.4) do not cancel. It is easy to understand from the form factors axioms'^ that the fields cr^{x) are 
semilocal with respect to a^{x). There are two semilocal fields in the critical sine-Gordon theory e^"^/^ 

with 'fi{x) = jUac '^y^ ^ofiy) ■ So we can suppose that 

ct±(.t) - e±*'^(^) . 



5. Discussion 

We conclude that the explicit form of the form factors in the scaling limit of the six vertex model is consistent 
with Lukyanov's conjecture on the free field representation of the form factors of the exponential fields in 
the sine-Gordon model ioi- f3 — 1, a — 1/2. A plausible argumentation extends this consistency to arbitrary 
values of a. This investigation can be continued in two ways. First, it would be interesting to calculate 
form factors of more general local fields, e. g. products of adjacent polarizations. Second, it is important to 
generalize this derivation to other models. For example, the bosonization of the Andrews-Baxter-Forrester 
model must give in the scaling limit the bosonization for the restricted sine-Gordon model. 

We expect that the divergency in the integral in /x is a peculiarity of the six vertex model or of the 
critical sine-Gordon model. It is related to the logarithmic factors in the scaling behavior of quantities, and 
originates in the dynamical mass generation phenomenon. 
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